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We study the Metal-Insulator transition in one-dimensional Hubbard superlattices (SL's), modelled 
by a repeated pattern of repulsive (i.e., positive on-site coupling) and free sites. The evolution of 
the local moment and of the charge gap (calculated from Lanczos diagonalization of chains up to 18 
sites), together with a strong coupling analysis, show that the electron density at which the system 
is insulating increases with the size of the free layer, relative to the repulsive one. In the insulating 
state, the mechanism of interaction between fermions separated by a free layer is the analog of 
superexchange, and the charge gap display universal features. 
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After many decades, the Metal-Insulator transition 
(MIT) still poses a great challenge, from both theoret- 
ical and experimental points of view ffl. For instance, 
one of the fundamental difficulties with MIT's is the 
identification of an order parameter which, in analogy 
with other (e.g., magnetic) phase transitions, would lead 
to spatial- dependent correlation functions and correla- 
tion lengths 0)- Nevertheless, considerable progress has 
been achieved in understanding and describing this phe- 
nomenon. MIT's can be driven by electronic correlations 
(the Mott transition), by disorder (Anderson transition), 
or by structural instabilities. The Mott transition, which 
is the one we discuss here, can take place as one varies 
any of the thermodynamic parameters, such as electron 
density, temperature, or the strength of correlation en- 
ergy. 

In the quest for a deeper understanding, it is impor- 
tant to gather as much information as possible about 
systems undergoing a MIT, by means of simplified mod- 
els amenable to well-controlled approximations, and the 
Hubbard model provides a nice testing ground for new 
ideas. In effect, the ground state of the half-filled homo- 
geneous chain corresponds to an insulating Spin-Density- 
Wave (SDW) state, which resembles the Neel state, 
but with power-law decay of correlations ||; its strong- 
coupling limit describes localized spins coupled antiferro- 
magnetically through a Heisenberg exchange interaction. 
Away from half-filling the ground state is metallic, with 
weak SDW correlations In two dimensions, the con- 
current van Hove singularity and nesting of the Fermi sur- 
face at half-filling give rise to an insulating antiferromag- 
netic (AFM) ground state for any finite on-site Coulomb 
repulsion U . As one dopes away from half-filling, the 
system undergoes a transition to a paramagnetic metal 
|^],[|, with weak short-ranged incommensurate spin-spin 
correlations |t]|| . Modifying the band structure through 
the addition of a second-neighbor hopping term does not 
seem to displace the transition point from half- filling H . 



Considerable insight into the MIT problem should there- 
fore be gained by seeking instances in which such dis- 
placement can occur. 

Motivated by the oscillatory exchange coupling in mag- 
netic metallic multilayers pj, we have recently stud- 
ied the effects of electronic correlations in superlattices 
(SL's), through a one-dimensional Hubbard- like model. 
In spite of its simplicity, a number of remarkable fea- 
tures were found, in marked contrast with the otherwise 
homogeneous system jll]] : Local moment weight could be 
transferred from repulsive to free sites; SDW quasi-order 
was wiped out as a result of frustration; and strong SDW 
correlations (in a subset of sites) could set in above half- 
filling. Since the combination of strong correlations with 
SL structures can modify the magnetic properties in such 
non-trivial way, transport properties should, most likely, 
be also affected. With this in mind, here we investigate 
the MIT transition in this one-dimensional Hubbard su- 
perlattice model. 

The model consists of a periodic arrangement of Ljj 
sites ("layers") in which the on-site coupling is repulsive, 
followed by L$ free {i.e., U — 0) sites, such that the 
Hamiltonian is written as 

U = -t ^2 ( c i<T c m<T + H.c.J + ^2 U * n A n H C 1 ) 

i, a i 

where, in standard notation, i runs over the sites of a 

t 

one-dimensional lattice, c ia (ci a ) creates (annihilates) a 

fcrmion at site i in the spin state a =f or |, and 

I 

n; = n,f + riii, w ith ni a = c ia .Ci a . The on-site Coulomb 
repulsion is taken to be site-dependent: Ui = U > 0, 
for sites within the repulsive layers, and Ui — other- 
wise; throughout this paper, numerical values of U will 
be given in units of t. It is important to notice that the SL 
structure breaks particle-hole symmetry even on bipar- 
tite lattices, as it can be seen through a simple Hartree- 
Fock argument: The exchange splitting is zero on free 
sites and non-zero on the repulsive ones, thus giving rise 



1 



to a non-uniform shift in the symmetric one-body local 
density of states. Differently from the homogeneous sys- 
tem, particle-hole symmetry cannot be restored by a uni- 
form shift in the chemical potential. 

The MIT in this model is studied numerically in the 
canonical ensemble. In addition to the numbers of sites, 
N s , and electrons, N e , another important parameter is 
the number of periodic cells, N c = N s /Nb, for a ba- 
sis with Nb — Ljj + Lq sites; in some cases, we were 
able to reach lattices with N s = 18 sites. The ground 
state \ipo) an d energy of Eq. (R) are obtained with the 
aid of the Lanczos algorithm ]12|-|l4[]: Starting with a 
trial state, the Hamiltonian is used to generate a second 
state, orthogonal to the first, so one ends up with a 2 x 2 
representation for the Hamiltonian ]ia). The diagonal- 
ization is trivial, leading to an estimate for the ground 
state and energy, which are used as new inputs for the 
subsequent iteration. This process is repeated until nu- 
merical convergence for the ground state energy has been 
achieved. In order to improve extrapolations to the ther- 
modynamic limit, boundary conditions were chosen such 
that the Fermi momentum was always one of the allowed 
k- values . The influence on the results can be assessed 
by simultaneously fixing system size, occupation and SL 
configuration, and performing calculations for different 
boundary conditions: the ground state energy turned out 
to be quite insensitive to the condition imposed (at most 
of the order of 1%) and similarly for the local moment. 
When dealing with the charge gap (see below) , the above 
mentioned choice has the additional advantage of provid- 
ing a trend (with system size) free from oscillations [Gl| . 

We assess whether the system is metallic or insulator 
by calculating the local moment and the charge gap, as 
we now discuss in turn. The local moment at site i is 
defined as 



(2) 



where m, = — n^; ensemble averages should be un- 
derstood as ground state averages, (...) = (ipo\ • ■ ■ \ipo)- 
Being a measure of both magnetism and degree of itin- 
erancy of the system, the local moment is useful in the 
investigation of the metallic or insulating character of 
the ground state. Indeed, in the case of a homogeneous 
lattice, for a finite value of the on-site repulsion U, the 
local moment is sharply peaked at half-filling; in the com- 
pletely localized limit, U = oo, (S 2 ) = 3/4. Above half 
filling, (>Sf ) decreases due to an increase in the double oc- 
cupancy of the sites. Unlike the magnetization (or even 
the sublattice magnetization), (m^), which vanishes iden- 
tically on a finite system due to the lack of spontaneous 
symmetry breaking, the local moment is always non-zero 
on a homogeneous lattice, except for p = and 2. 

For a SL the site-dependent Coulomb repulsion leads 
to a non-uniform distribution of local moments through- 
out the lattice Oj . Figure [l] shows typical plots of the 



density dependence of the local moment at a repulsive 
site, for lattices with N s = 12. As the SL configuration 
is changed, the maximum value of (Sf) still approaches 
3/4 as U is increased; its position is displaced to higher 
fillings, without showing any [/-dependence. By analogy 
with the homogeneous case, one is led to identify this 
peak position with the density, pi, at which the system 
becomes an insulator; Fig. [jjthen shows that pj increases 
continuously with Lq, for fixed Ljj . To see how this comes 
about, we should ask ourselves how can the SL be filled 
up with electrons in a way to obtain maximum hopping 
hindering and largest moment on the repulsive sites. In 
strong coupling this is equivalent to pin the electrons, 
and is achieved by placing two on each of the free sites 
(thus rendering them magnetically inert) and one on each 
of the repulsive sites (maximum polarization); this leads 
to 



Pi = 



2L + L, 



A J 



(3) 



The maxima position of Fig. [l] are given exactly by Eq. 
(^), with Ljj = 1, and Lq = 1, 2, 3, and 5. Configurations 
with Ljj > 1 follow the same pattern, and reflect the fact 
that the layer thicknesses come into the definition of pi 
only through the ratio Ljj/Lq. Also, note that in the 
limit of a uniform lattice, Lq — > 0, one has pi — > 1, thus 
recovering the insulating behavior exactly at half filling. 
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FIG. 1. Local moment at a repulsive site versus p for 12 
site chains with Ljj — 1 and Lo = 1 (a), 2 (b), 3 (c), and 
5 (d). Squares and dotted lines respectively represent data 
for (7 = and oo for the homogeneous system; triangles and 
circles respectively represent data for superlattices with (7 = 3 
and 12. 



The above discussion allows us to perform a strong cou- 
pling perturbation theory analysis |l6| . Let us start with 
the dimerized SL configuration; the relevant Hilbert sub- 
space is the one with repulsive and free sites being respec- 
tively singly and doubly occupied. The first non-trivial 
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contribution from the hopping term comes in fourth or- 
der, and gives rise to an effective Heisenberg Hamiltonian 



[k,e] 



(4) 



where J2 = 8t 4 /\U\ 3 , and k, £ are repulsive sites (i.e., 
next-nearest neighbors on the original dimerized lat- 
tice); trivial constants have been dropped from Eq. (Q). 
Thus, the (magnetically inert) free sites intermediate the 
Heisenberg exchange interaction between moments on 
the repulsive sites; the analogy with the usual superex- 
change mechanism involving non-magnetic atoms arises 
here quite naturally. These arguments can be straight- 
forwardly generalized to other configurations, with suit- 
able changes: (i) within a repulsive layer (Ljj > 1), 
spins interact via the usual second order Heisenberg ex- 
change coupling Ji = 4t 2 /\U\; (ii) the first non-trivial 
coupling across a free layer of thickness Lq only ap- 
pears in 2(Lq + l)-th order, and the exchange coupling 
becomes ~ f 2 ( L o+ 1 ) /\U\ 2Lo+1 . The ensuing Heisenberg 
model is again defined only on the repulsive sites, with 
non-uniform, but periodic, exchange couplings. Away 
from strong coupling, numerical data are consistent with 
this picture 

Further evidence in favor of the MIT being located at 
pi comes from the analysis of the charge gap, defined by 

A c = E(N C , N e + l) + E(N C , N e -1)- 2E(N C , N e ), 

(5) 

where, for a given SL configuration, E(N c ,N e ) is the 
ground state energy for a chain with N c periodic cells 
and N e electrons. This gap is a measure of the charge 
excitation spectrum (i.e., the cost of adding extra par- 
ticles): When extrapolated to the thermodynamic limit, 
a non-zero value indicates an insulating state, whereas 
a zero value corresponds to a metallic state. Before dis- 
cussing our results, a comment is in order at this point. 
The Drude weight, which is proportional to the the sec- 
ond derivative of the ground state energy with respect 
to the magnetic flux through the ring [g7|-Q, is an el- 
egant test of the metallic or insulating character of the 
system: It yields impressive results when Bethe-Ansatz 
data (corresponding to arbitrarily large lattices) are used 
|l8|| . For small systems, however, extrapolations with the 
charge gap appear to be more conclusive than those with 
the Drude weight 

Figure |^(a) shows charge gap data for the dimerized 
configuration, Ljj = Lq = 1. In this case pi = 3/2, and 
extrapolation (through least-squares fits) of data towards 
N c — ► 00 yields finite values for A c , for both U = 6 and 
12, confirming the insulating character of this filling fac- 
tor. By contrast, for p — 5/3, the data extrapolate to 
A c ~ 0, for both values of U, signaling a metallic state; 
for clarity, the only metallic case shown in Fig. |2|(a) cor- 
responds to p = 5/3, but this is a representative example 



for all p 7^ pi. For the configuration Ljj = 1, Lq = 2, pi 
is now 5/3, and Fig. ||(b) shows that the roles of fillings 
5/3 and 3/2 have been reversed with respect to that of 
Fig. |(a). 
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FIG. 2. Charge gap vs. the inverse number of cells for the 
SL configurations Lu = Lo = 1 (a) and Lu — 1, Lo = 2 (b), 
and densities p — 3/2 (empty symbols) and p = 5/3 (filled 
symbols). Circles and squares denote U = 6 and U = 12, 
respectively. The lines are least-squares fits to the data. 
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FIG. 3. Scaled extrapolated-charge-gap versus scaled 
on-site energy. Data for the configurations (Lu,Lo) are rep- 
resented as follows: (1,1) by squares; (1,2) by circles; (1,3) by 
up-triangles; and (2,2) by down-triangles. The full line is the 
Mott-Hubbard gap for the homogeneous system, calculated 
from the Lieb-Wu solution. The inset shows the intermediate 
coupling region. 



Similarly to the homogeneous system, the gap in- 
creases with U. This similarity, however, goes beyond a 
simple qualitative trend: If both the extrapolated charge 
gap and the on-site coupling U are scaled by pj, strong 
coupling data for different pi fall on a universal function, 
as shown in Fig. ^|. This function, in turn, reduces to the 
Mott-Hubbard gap for the homogeneous system (pi = 1), 
known for any U through the Bethe ansatz solution [p|,p2| ; 
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its asymptotic forms are A c ~ (8\/tU /n) exp (—2-irt/U), 
in weak coupling, and A c ~ U, in strong coupling. As 
U decreases, the relative errors in extrapolating the gaps 
increase; nonetheless, the resulting fits to the scaled Lieb- 
Wu gap (shown as the inset of Fig. ||) are still satisfied 
within error bars. 

In summary, within a simple model wc have estab- 
lished that the superlattice structure induces a shift in 
the critical density at which the Mott-Hubbard transi- 
tion occurs: It increases with the ratio between free and 
repulsive layer thicknesses. This comes about as a result 
of having to doubly (singly) occupy the free (repulsive) 
sites, in order to keep electronic motion to a minimum, 
i.e., to that resulting from quantum fluctuations. These 
insulating states have interesting properties: (i) fermions 
on free sites act as mediators of the magnetic exchange 
interaction between fermions on repulsive sites ('superex- 
change'); (ii) if properly scaled, the charge gap seems to 
follow a universal curve (reducing to that for the homoge- 
neous system at half filling) , for all superlattices. Exper- 
imental realizations of (effectively) one-dimensional sys- 
tems undergoing Mott-Hubbard transitions are, by now, 
well established in organic compounds p!| . In Bcchgaard 
salts, for instance, organic complexes are stacked along 
a given direction; their large anisotropy is responsible 
for the one-dimensional character of the system. Within 
closely related families of these compounds, some are in- 
sulators whereas others are metals, and an interesting 
possibility would be to stack layers of these different com- 
plexes, forming a one-dimensional superlattice, and in- 
vestigate the evolution of the metal-insulator transitions 
with the superlattice structure; we hope the scenario pro- 
posed here stimulates experimental studies along these 
lines. 
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